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The properties of cosmic rays with energies above 10 GeV have to be deduced 
from the spacetime structure and particle content of the air showers which they 
initiate. In this review, a summary of the phenomenology of these giant air showers 
is presented. We describe the hadronic interaction models used to extrapolate 
results from collider data to ultra high energies, an also the main electromagnetic 
processes that govern the longitudinal shower evolution as well as the lateral spread 
of particles. 



1. Introduction 

For primary cosmic ray energies above 10® GeV, the flux becomes so low 
that direct detection of the primary particle using detectors in or above 
the upper atmosphere is not longer possible. In these cases the primary 
particle has enough energy to initiate an extensive air shower (EAS) in 
the atmosphere. If the primary cosmic ray particle is a nucleon or nucleus 
the cascade begins with a hadronic interaction. The number of hadrons 
increases through subsequent generations of particle interactions. In each 
generation about 20% of the energy is transferred to an electromagnetic 
cascade by rapid decays of neutral pions. Ultimately, the electromagnetic 
cascade dissipates roughly 90% of the primary particle's energy trough ion- 
ization. The remaining energy is carried by muons and neutrinos from 
charged pion decays. The electromagnetic and weak interactions are rather 
well understood. However, uncertainties in hadronic interactions at ultra 
high energies constitute one of the most problematic sources of systematic 
error in analysis of air showers. In what follows a brief report of the phe- 
nomenology of EAS is presented, focusing in those aspects which are the 
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main source of systematic uncertainties affecting somehow the determina- 
tion of primary energy and mass composition. A complete review on the 
phenomenology of cosmic ray air showers can be found in ^. 

2. Hadronic Processes 

Soft nniltiparticle production with small transverse momenta with respect 
to the collision axis is a dominant feature of most hadronic events at center- 
of-mass energies 10 GeV < ^/s < 50 GeV. Despite the fact that strict 
calculations based on ordinary QCD perturbation theory are not feasible, 
there are some phenomenological models that successfully take into ac- 
count the main properties of the soft diffractive processes. These models, 
inspired by 1/N QCD expansion are also supplemented with generally ac- 
cepted theoretical principles like duality, unitarity, Regge behavior, and 
parton structure. The interactions arc no longer described by single parti- 
cle exchange, but by highly complicated modes known as Reggeons. Up to 
about 50 GeV, the slow growth of the cross section with ^/s is driven by a 
dominant contribution of a special Reggcon, the Pomcron. 

At higher energies, semihard interactions arising from the hard scatter- 
ing of partons that carry only a very small fraction of the momenta of their 
parent hadrons can compete successfully with soft processes. These semi- 
hard interactions lead to the minijet phenomenon, i.e., jets with transverse 
energy {Et = \Pt\) much smaller than the total center-of-mass energy. Un- 
like soft processes, this low-p^ jet physics can be computed in perturbative 
QCD. The parton-parton minijet cross section is given by 

aQCD(s,pr°") = E / ^ d\i\^X,Mx^,\i\) X,Mx,,\i\) 

Xl J X2 JqI^.^ d\t\ 

(1) 

where xi and X2 are the fractions of the momenta of the parent hadrons 
carried by the partons which collide, d(Jij/d\t\ is the cross section for 
scattering of partons of types i and j according to elementary QCD di- 
agrams, fi and fj are parton distribution functions (pdf's), s = xix^s 
and —t = s(l — cos'i?*)/2 = are the Mandelstam variables for this 
parton-parton process, and the sum is over all parton species. The inte- 
gration limits satisfy Q^j^ < \t\ < s/2, with Qmin the minimal momentum 
transfer. 

A first source of uncertainty in modeling cosmic ray interactions at 
ultra high energy is encoded in the extrapolation of the measured parton 
densities several orders of magnitude down to low x. Primary protons that 
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impact on the upper atmosphere with energy w 10^^ GeV, yield partons 
with X = 2p*/^ « rnT^ I ^fa ^ 10^^, whereas current data on quark and 
gluon densities are only available for x ~ 10~^ to within an experimental 
accuracy of 3% for ~ 20 GeV^ ^. Moreover, application of HERA data 
to baryonic cosmic rays assumes universality of the pdf's. 

For large and not too small x, the Dokshitzcr-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations successfully predict the dependence 
of the quark and gluon densities (g' and g, respectively). In the double- 
leading-logarithmic approximation the DGLAP equations predict a steeply 
rising gluon density, xg ~ x~^^'^ ^ which dominates the quark density at low 
X, in agreement with experimental results obtained with the HERA col- 
lider ^. Specifically, HERA data are found to bo consistent with a power 
law, xg{x, Q^) ^ x~^^ ^ with an exponent Ah between 0.3 and 0.4 ^ . 

The high energy minijet cross section is then determined by the small-x 
behavior of the parton distributions or, rather, by that of the dominant 
gluon distribution (via the lower limits of the x\, X2 integrations) which 
gives ^: 

c7qcd(s)oc r ^xr^« r ^X2-^«~s^«ln(s/so), (2) 

where sq is a normalization constant. One caveat is that the inclusive 
QCD cross section given in Eq. (2) is a Born approximation, and therefore 
automatically violates unitarity. 

The procedure of calculating the inelastic cross section from inclusive 
cross sections is known as unitarization. In the eikonal model of high 
energy hadron-hadron scattering,the inelastic cross section, assuming a real 
eikonal function, is given by 

cTinei = j(fb[l-exp [-2x,,« {s,b)- 2x,,,, (s, b)] } , (3) 

where the scattering is compounded as a sum of QCD ladders via hard and 
soft processes through the eikonals Xh^id ^^'^ Xaott • should be noted that 
we have ignored spin-dependent effects and the small real part of the scat- 
tering amplitude, both good approximations at high energies. Now, if the 
eikonal function, x(.s, b) = Xaon (*i ^) + Xhard i-^i ^) ~ -^/2, indicates the mean 
number of partonic interaction pairs at impact parameter b, the probability 
Pn for having n independent partonic interactions using Poisson statistics 
reads, Pn = (A"/n!)e~'^. Therefore, the factor 1 — e^^^ = J2'^=iPn in 
Eq. (3) can be interpreted scmiclassically as the probability that at least 1 
of the 2 protons is broken up in a collision at impact parameter b. With 
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Figure 1. Energy dependence of the pp inelastic cross section as predicted by Eqs. (5) 
and (6) with 0.3 < Ajj < 0.4. The darkly shaded region between the solid lines corre- 
sponds to the model with Gaussian parton distribution in b. The region between the 
dashed-dotted lines corresponds to the model with exponential fall-off of the parton 
density in b. 

this in mind, the inelastic cross section is simply the integral over all colli- 
sion impact parameters of the probability of having at least 1 interaction, 
yielding a mean minijet multiplicity of (njet) ~ crqcD/cinci The lead- 
ing contenders to approximate the (unknown) cross sections at cosmic ray 
energies, SIBYLL and QGSJET share the eikonal approximation but 
differ in their ansdtse for the eikonals. In both cases, the core of dominant 
scattering at very high energies is the parton-parton minijet cross section 
given in Eq. (1), 

X.„a = ^'^QCD(s,Pr°^) (4) 

where the normalized profile function, j d^b A(s,b) = 1, indicates the dis- 
tribution of partons in the plane transverse to the collision axis. 

In the QGSJET-like models, the core of the hard eikonal is dressed with 
a soft-pomeron pre-evolution factor. This amounts to taking a parton dis- 
tribution which is Gaussian in the transverse coordinate distance |&|. 

In SiBYLL-like models, the transverse density distribution is taken as 
the Fourier transform of the proton electric form factor, resulting in an 
energy-independent exponential (rather than Gaussian) fall-off of the par- 



February 2, 2008 18:31 Proceedings Trim Size: 9in x 6in 



procfaro 



5 




Figure 2. Left panel: The slowly rising curves indicate the mean inelasticity in proton 
air collisions as predicted by QGSJET and SIBYLL. The falling curves indicate the proton 
mean free path in the atmosphere. Right panel: Mean multiplicity of charged secondary 
particles produced in inelastic proton-air collisions processed with QGSJET and SIBYLL. 

ton density profile with \b\. The main characteristics of the pp cascade 
spectrum resulting from these choices are readily predictable: the harder 
form of the SIBYLL form factor allows a greater retention of energy by the 
leading particle, and hence less available for the ensuing shower. Conse- 
quently, on average SiBYLL-like models predict a smaller multiplicity than 
QGSJET-like models (see e.g. I6,i7,i8,i9-)^ 

At high energy, x^ast ^ Xhard' inelastic cross section is dom- 

inated by the hard eikonal. With the appropriate choice of normalization, 
the cross section in Eq. (2) can be well-approximated by a power law. This 
implies that the growth of the inelastic cross section according to QGSJET- 
like models is given by 

CTinci J d% Q{hs-\h\) = ^hl ^ 4^a^ff Ah In'ls/so) 0.52 Ah \n{s/so) mb . 

(5) 

For SiBYLL-like models, the growth of the inelastic cross section also satu- 
rates the In^ s Froissart bound, but with a multiplicative constant which is 
larger than the one in QGSJET-like models Namely, 

CTinci 3.2 AHln^(s/so) mb . (6) 

Figure 1 illustrates the large range of predictions for pp inelastic cross 
section which remain consistent with HERA data. When the two leading 
order approximations discussed above are extrapolated to higher energies. 
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both are consistent with existing cosmic ray data. Note, however, that 
in both cases the range of allowed cross-sections at high energy varies by 
a factor of about 2 to 3. The points in Fig. 1 correspond to the most 
up-to-date estimate of the pp cross section from cosmic ray shower data . 

There are three event generators, SIBYLL QGSJET and DPMJET 
which are tailored specifically for simulation of hadronic interactions up to 
the highest cosmic ray energies. The latest versions of these packages are 
SIBYLL 2.1 QGSJET 01 and DPMJET III "^"^ ; respectively. In QGSJET, 
both the soft and hard processes are formulated in terms of Pomeron ex- 
changes. To describe the minijets, the soft Pomeron mutates into a "semi- 
hard Pomeron" , an ordinary soft Pomeron with the middle piece replaced 
by a QCD parton ladder, as sketched in the previous paragraph. This is 
generally referred to as the "quasi-eikonal" model. In contrast, SIBYLL and 
DPMJET follow a "two channel" eikonal model, where the soft and the semi- 
hard regimes are demarcated by a sharp cut in the transverse momentum: 
SIBYLL 2.1 uses a cutoff parametrization inspired in the double leading log- 
arithmic approximation of the DGLAP equations, whereas dpmjet III uses 
an ad hoc parametrization for the transverse momentum cutoff ^. 

The transition process from asymptotically free partons to colour- 
neutral hadrons is described in all codes by string fragmentation models 
Different choices of fragmentation functions can lead to some differences in 
the hadron multiplicities. However, the main difference in the predictions 
of QGSJET-like and SiBYLL-like models arises from different assumptions in 
extrapolation of the parton distribution function to low energy. 

Now we turn to nucleus-nucleus interactions, which cause additional 
headaches for event generators which must somehow extrapolate pp inter- 
actions in order to simulate the proton- air collisions of interest. All the 
event generators described above adopt the Glauber formalism 

Since the codes described above are still being refined, the disparity 
between them can vary even from version to version. At the end of the 
day, however, the relevant parameters boil down to two: the mean free 
path, A = (cTprod ?^)~^, and the inelasticity, K = 1 — Eiead/Eproj, where n 
is the number density of atmospheric target nucleons, -Eioad is the energy 
of the most energetic hadron with a long lifetime, and -Bproj is the energy 
of the projectile particle. Overall, SIBYLL has a shorter mean free path 
and a smaller inelasticity than qgsjet, as indicated in Fig. 2. Since a 
shorter mean free path tends to compensate a smaller inelasticity, the two 
codes generate similar predictions for an air shower which has lived through 
several generations. The different predictions for the mean charged particle 



February 2, 2008 18:31 Proceedings Trim Size: 9in x 6in 



procfaro 



7 

multiplicity in proton-air collisions are shown in Fig. 2. Both models predict 
the same multiplicity below about 10'' GeV, but the predictions diverge 
above that energy. Such a divergence readily increases with rising energy. 
As it is extremely difficult to observe the first interactions experimentally, 
it is not straightforward to determine which model is closer to reality. 

3. Electromagnetic Component 

The evolution of an extensive air shower is dominated by electromagnetic 
processes. The interaction of a baryonic cosmic ray with an air nucleus high 
in the atmosphere leads to a cascade of secondary mesons and nucleons. The 
first few generations of charged pions interact again, producing a hadronic 
core, which continues to feed the electromagnetic and muonic components of 
the showers. Up to about 50 km above sea level, the density of atmospheric 
target nucleons is n ~ 10^'^ cm~^, and so even for relatively low energies, 
say E^± « 1 TeV, the probability of decay before interaction falls below 
10%. Ultimately, the electromagnetic cascade dissipates around 90% of the 
primary particle's energy, and hence the total number of electromagnetic 
particles is very nearly proportional to the shower energy. 

By the time a vertically incident 10"'^'^ GcV proton shower reaches the 
ground, there are about 10^^ secondaries with energy above 90 keV in the 
the annular region extending 8 m to 8 km from the shower core. Of these, 
99% are photons, electrons, and positrons, with a typical ratio of 7 to e+e" 
of 9 to 1. Their mean energy is around 10 MeV and they transport 85% 
of the total energy at ground level. Of course, photon-induced showers are 
even more dominated by the electromagnetic channel, as the only significant 
muon generation mechanism in this case is the decay of charged pions and 
kaons produced in 7-air interactions. 

It is worth mentioning that these figures dramatically change for the case 
of very inclined showers. For a primary zenith angle, 9 > 70°, the electro- 
magnetic component becomes attenuated exponentially with atmospheric 
depth, being almost completely absorbed at ground level. We remind the 
reader that the vertical atmosphere is « 1000 g/cm^, and is about 36 times 
deeper for completely horizontal showers. 

In contrast to hadronic collisions, the electromagnetic interactions of 
shower particles can be calculated very accurately from quantum electro- 
dynamics. Electromagnetic interactions are thus not a major source of 
systematic errors in shower simulations. The first comprehensive treat- 
ment of electromagnetic showers was elaborated by Rossi and Greissen 
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This treatment was recently cast in a more pedagogical form by Gaisser 
and a summary is presented in ^ . 

The generation of the electromagnetic component is driven by electron 
bremsstrahlung and pair production . Eventually the average energy per 
particle drops below a critical cncrgV: cq, at which point ionization takes 
over from bremsstrahlung and pair production as the dominant energy loss 
mechanism. The energy loss rate due to bremsstrahlung radiation is 
nearly proportional to their energy, whereas the ionization loss rate varies 
only logarithmically with the energy. The changeover from radiation 
losses to ionization losses depopulates the shower. One can thus categorize 
the shower development in three phases: the growth phase, in which all the 
particles have energy > eq; the shower maximum, Xmaxi and the shower 
tail, where the particles only lose energy, get absorbed or decay. 

The relevant quantities participating in the development of the electro- 
magnetic cascade are the probability for an electron of energy E to radiate 
a photon of energy k = yE and the probability for a photon to produce a 
pair e+e~ in which one of the particles (hereafter e~) has energy E = xk. 
These probabilities are determined by the properties of the air and the cross 
sections of the two processes. 

In the energy range of interest, the impact parameter of the electron or 
photon is larger than an atomic radius, so the nuclear field is screened by 
its electron cloud. In the case of complete screening, where the momentum 
transfer is small, the cross section for bremsstrahlung can be approximated 



where Af.s is the effective mass number of the air, Xq is a constant, and 
Na is Avogadro's number. In the infrared limit (i.e., y <C 1) this approxi- 
mation is inaccurate at the level of about 2.5%, which is small compared to 
typical experimental errors associated with cosmic air shower detectors. Of 
course, the approximation fails as y ^ 1, when nuclear screening becomes 
incomplete, and as y — > 0, at which point the LPM and dielectric suppres- 
sion effects become important. This infrared divergence is eliminated by the 
interference of bremsstrahlung amplitudes from multiple scattering centers. 
This collective effect of the electric potential of several atoms is known as 
the Landau-Pomeranchuk-Migdal (LPM) effect 33,34 Using similar approx- 
imations, the cross section for pair production can be obtained . 

The LPM suppression of the cross section results in an effective increase 
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Figure 3. Left panehAverage longitudinal shower developments of 10^^ GeV proton 
(dashed-dotted line) and 7-rays with and without the LPM effect (solid and dotted lines, 
respectively). The primary zenith angle was set to Q = 60°. Right panel: Longitudinal 
development of muons and electrons as a function of the slant depth for lO'^^ GeV 
proton-induced showers. 

of the mean free path of electrons and photons. This effectively retards the 
development of the electromagnetic component of the shower. It is natural 
to introduce an energy scale, £'lpm, at which the inelasticity is low enough 
that the LPM effect becomes significant 

The experimental confirmation of the LPM effect at Stanford Linear 
Accelerator Center (SLAC) has motivated new analyses of its conse- 
quences in cosmic ray physics 37,38,39,40 _ most evident signatures of 
the LPM effect on shower development are a shift in the position of the 
shower maximum Xmax and larger fluctuations in the shower development. 

Since the upper atmosphere is very thin the LPM effect becomes notice- 
able only for photons and electrons with energies above i?LPM ~ 10^" GeV. 
For baryonic primaries the LPM effect does not become important until the 
primary energy exceeds lO^^GeV. To give a visual impression of how the 
LPM effect slows down the initial growth of high energy photon-induced 
showers, we show the average longitudinal shower development of I0^° GeV 
proton and 7-ray showers (generated using aires 2.6.0 '*^) with and without 
the LPM effect in Fig. 3. 

At energies at which the LPM effect is important (wzz., E > i?LPM), 
7-ray showers will have already commenced in the geomagnetic field at al- 
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most all latitudes. This reduces the energies of the primaries that reaeh the 
atmosphere, and thereby compensates the tendency of the LPM effect to 
retard the shower development. The first description of photon interactions 
in the geomagnetic field dates back at least as far as 1966 with a punctu- 
ated revival of activity in the early 1980's More recently, a rekindling of 
interest in the topic has led to refined calculations ^4,45 pj-jniary photons 
with energies above 10^° GeV convert into e+e~ pairs, which in turn emit 
synchrotron photons. Regardless of the primary energy, the spcctnim of 
the resulting photon "prcshower" entering the upper atmosphere extends 
over several decades below the primary photon energy, and is peaked at 
energies below 10^'^ GeV The geomagnetic cooling thus switches on at 
about the same energy at which the LPM effect does, and thereby preempts 
the LPM-related observables which would otherwise be evident. 

The relevant parameter to determine both conversion probability and 
synchrotron emission in E x B± , where E is the 7-ray energy and B± the 
transverse magnetic field. This leads to a large directional and geographi- 
cal dependence of shower observables. Thus, each experiment has its own 
preferred direction for identifying primary gamma rays. 

3.0.1. Electron lateral distribution function 

The transverse development of electromagnetic showers is dominated by 
Coulomb scattering of charged particles off the nuclei in the atmosphere. 
The lateral development in electromagnetic cascades in different materials 
scales well with the Moliere radius Tm = Eg Xq/eq, which varies inversely 

with the density of the medium, rM = rM{hoh) ^'^p'l'''\h) ^ — ^'p ^^(h) ' where 
Eg !v 21 MeV and the subscript OL indicates a quantity taken at a given 
observation level. 

Approximate calculations of cascade equations in three dimensions to 
derive the lateral structure function for a pure electromagnetic cascade in 
vertical showers were obtained by Nishimura and Kamata and later 
worked out by Greisen in the well-known NKG formula, 

pir)^^c(^) 1 + - , (8) 

where is the total number of electrons, r is the distance from the 
shower axis. For a primary of energy Eq, the so-called "age parameter", 
^NKG = <^/(l + ^ in(-Eo/eo) -j^ characterizes the stage of the shower devel- 
opment in terms of the depth of the shower in radiation lengths, i.e., 

t = Pa.tm{z) dz/Xo. 
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The NKG formula may also be extended to describe showers initiated 
by baryons In such an extension, one finds a deviation of behavior of 
the Moliere radius when using a value of the age parameter which is de- 
rived from theoretical predictions for pure electromagnetic cascades. It is 
possible to generalize the NKG formula for the electromagnetic component 
of baryon- induced showers by modifying the exponents in Eq. (8) The 
derived NKG formula provides a good description of the e+e~ lateral dis- 
tribution at all stages of shower development for values of r sufficiently far 
from the hadronic core. Fortunately, this is the experimentally interesting 
region, since typical ground arrays can only measure densities at r > 100 m 
from the shower axis, where detectors are not saturated. It should be 
mentioned that an NKG-like formula can be used to parametrize the total 
particle's density observed in baryon-induced showers 

In the case of inclined showers, one normally analyzes particle densities 
in the plane perpendicular to the shower axis. Simply projecting distribu- 
tions measured at the ground into this plane is a reasonable approach for 
near- vertical showers, but is not sufficient for inclined showers. In the latter 
case, additionally asymmetry is introduced because of both unequal atten- 
uation of the electromagnetic components arriving at the ground earlier 
than and later than the core Moreover, deflections on the geomagnetic 
field become important for showers inclined by more than about 70°. 

In the framework of cascade theory, any effect coming from the influence 
of the atmosphere should be accounted as a function of the slant depth 
t Following this idea, a LDF valid at all zenith angles < 70° can be 
determined by considering 

t'ieX) =t sec0il + K coaCy\ (9) 

where ( is the azimuthal angle in the shower plane, K = Kq ta,n9, and Kq 
is a constant extracted from the fit ^3,55 Tiign^ the particle lateral distri- 
butions for inclined showers p{r,t') are given by the corresponding vertical 
LDF p{r,t) but evaluated at slant depth t'{9,Q where the dependence on 
the azimuthal angle is evident. 

For zenith angles > 70°, the surviving electromagnetic component 
at ground is mainly due to muon decay and, to a much smaller extent, 
hadronic interactions, pair production and bremsstrahlung. As a result 
the lateral distribution follows that of the muon rather closely. In Fig. 3 
the longitudinal development of the muon and electron components are 
shown. It is evident from the figure that for very inclined showers the 
electromagnetic development is due mostly to muon decay 
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3.1. The muon component 

The muonic component of EAS differs from the electromagnetic compo- 
nent for two main reasons. First, muons are generated through the decay 
of cooled charged pions, and thus the muon content is sensitive to the ini- 
tial baryonic content of the primary particle. Furthermore, since there is 
no "muonic cascade" , the number of muons reaching the ground is much 
smaller than the number of electrons. Specifically, there are about 5 x 10® 
muons above 10 MeV at ground level for a vertical 10^^ GeV proton induced 
shower. Second, the muon has a much smaller cross section for radiation 
and pair production than the electron, and so thc^ miionic component of 
EAS develops differently than does the electromagnetic component. The 
smaller multiple scattering suffered by muons leads to earlier arrival times 
at the ground for muons than for the electromagnetic component. 

The ratio of electrons to muons depends strongly on the distance from 
the core; for example, the e+e~ to iX^n~ ratio for a 10^^ GeV vertical proton 
shower varies from 17 to 1 at 200 m from the core to 1 to 1 at 2000 m. The 
ratio between the electromagnetic and muonic shower components behaves 
somewhat differently in the case of inclined showers. For zenith angles 
greater than 60°, the e+e~//i+/u~ ratio remains roughly constant at a given 
distance from the core. As the zenith angle grows beyond 60°, this ratio 
decreases, until &t 9 = 75°, it is 400 times smaller than for a vertical 
shower. Another difference between inclined and vertical showers is that 
the average muon energy at ground changes dramatically. For horizontal 
showers, the lower energy muons are filtered out by a combination of energy 
loss mechanisms and the finite muon lifetime: for vertical showers, the 
average muon energy is 1 GeV, while for horizontal showers it is about 2 
orders of magnitude greater. 

High energy muons lose energy through e+e~ pair production, muon- 
nucleus interaction, bremsstrahlung, and knock-on electron ((5-ray) pro- 
duction The first three processes are discrete in the sense that they 
are characterized by high inelasticity and a large mean free path. On the 
other hand, because of its short mean free path and its small inelasticity, 
knock-on electron production can be considered a continuous process. The 
muon bremsstrahlung cross section is suppressed by a factor of (we/m^)^ 
with respect to electron bremsstrahlung, see Eq. (7). Since the radiation 
length for air is about 36.7 g/cm^, and the vertical atmospheric depth is 
1000 g/cw?, muon bremsstrahlung is of negligible importance for vertical 
air shower development. Energy loss due to muon- nucleus interactions is 
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somewhat smaller than muon brcmsstrahlung. Energy loss by pair produc- 
tion is slightly more important than bremsstrahlung at about 1 GeV, and 
becomes increasingly dominant with energy. Finally, knock-on electrons 
have a very small mean free path, but also a very small inelasticity, so that 
this contribution to the energy loss is comparable to that from the hard 
processes. 

In addition to muon production through charged pion decay, photons 
can directly generate muon pairs, or produce hadron pairs which in turn 
decay to muons. In the case of direct pair production, the large muon mass 
leads to a higher threshold for this process than for electron pair produc- 
tion. Furthermore, QED predicts that production is suppressed by a 
factor (TOe/m^)^ compared the Bcthc-Heitler cross section. The cross sec- 
tion for hadron production by photons is much less certain, since it involves 
the hadronic structure of the photon. This has been measured at HERA 
for photon energies corresponding to Ei-^h = 2 x 10* GeV. This energy is 
still well below the energies of the highest energy cosmic rays, but nonethe- 
less, these data do constrain the extrapolation of the cross sections to high 
energies. 

The muon content of the shower tail is quite sensitive to unknown de- 
tails of hadronic physics. This implies that attempts to extract composition 
information from measurements of muon content at ground level tend to 
be systematics dominated. The muon LDF is mostly determined by the 
distribution in phase space of the parent pions. However, the pionization 
process together with muon deflection in the geomagnetic field obscures the 
distribution of the first generation of pions. A combination of detailed sim- 
ulations, high statistics measurements of the muon LDF and identification 
of the primary species using uncorrelated observables could shed light on 
hadronic interaction models. 
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